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A quantum state contains the maximal amount of information available for a given quantum 
system. In this paper we use weak-value expressions to reconstruct quantum states of continuous- 
variable systems in the quantum optical domain. The role played by postselecting measured data 
will be particularly emphasized in the proposed setup, which is based on an interferometer just using 
simple homodyne detection. 
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I. INTRODUCTION 

Since the launch of the idea and formalism by 
Aharonov, Albert, and Vaidman 0, Q a considerable 
amount of research has focused on weak measurements. 
Basic questions behind the concept have been clarified in 
a series of papers [3j-|9| , and it was possible to extend the 
original proposal in various directions 

EMU- 

It is fasci- 
nating to see how the weak- value formalism touches very 
different fields like phase-space physi cs fTBI . [l6| . quan- 
tum trajectories [T7l.ll8j. contextuality [l9j|. and quantum 
cloning [2fj| . 

In parallel to these theoretical developments the ex- 
perimental realization pll ] of weak measurements has re- 
cently found a rich territory of applications [22I - I24I ] . Even 
the most fundamental debates about the realism of quan- 
tum observables [25I - t27l and the role of trajectories in 
double-slit experiments [28l. [29j can now be attacked with 
the assistance of weak measurements carried out in the 
laboratory. 

In the present paper we concentrate on yet another 
crosslink: Lately quantum state reconstruction [3(1 HU 
using weak measurements has become a topic of great 
interest [32 36]. Several works analyze the aspect of a 
weak interaction [37i - [39| in non-standard quantum state 
estimation. Obviously there is a special role played by 
postselection in creating weak values of observables. This 
thought is also pursued in Ref. 32] where a quantum 
wave function has been reconstructed using the special 
form of weak values. Later this procedure was general- 
ized for mixed quantum states [341 ]. drawing connections 
between the structure of weak values and certain phase- 
space distributions I40j- p2ll . A similar connection was also 
established in Refs. |43| and [lij and later used to build 
yet another reconstruction scheme [3(| [45[ . 

In this article we first recall the special form of weak 
values used to derive reconstruction relations. Our fo- 
cus then lies in establishing a direct connection between 
a quantum optical implementation and these reconstruc- 
tion relations. The main contribution of this work there- 
fore consists of working out a realizable quantum optical 



scheme that fully exploits the form of weak values to re- 
construct a continuous- variable state of a single mode of 
the quantized electromagnetic field. 

In Sec. |TT] we will recapitulate the formalism of weak 
measurements. Section IIIII provides the reconstruction 
relations for quantum states of continuous- variable sys- 
tems. In Sec. IIVI we show how to realize the method 
with a quantum optical setup. Finally, in Sec.|V]we sim- 
ulate the optical experiment for a special quantum state 
of light and discuss to which class of states the recon- 
struction scheme can be applied. 



II. WEAK MEASUREMENTS 

In this section we will briefly recall the concepts of 
weak measurements [l[ to introduce our notation and to 
make the paper self-consistent. An appropriate scaling 
allows us to choose dimensionless quantities everywhere. 

A weak measurement can be divided into four succes- 
sive steps as depicted in Fig. Q] In the first step the 
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FIG. 1. Basic setup of a weak measurement which can be 
applied to reconstruct the unknown state of a system S 
with the help of a meter M initially prepared in a specific 
state \(j>). After system and meter have been separately de- 
fined in the first step, they interact weakly in the second step. 
This leaves S and M in the joint state |$). After a suited pro- 
jective measurement with outcome bj in the third step, the 
postselected expectation value E' 3 '[om] of a convenient meter 
observable om can be obtained in the fourth step. We demon- 
strate in Sec. IIIII that the probability distribution of finding 
the eigenvalue bj and the expectation value E"' are sufficient 
to reconstruct the state lib). 
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system S is given in the unknown state \4>) and the me- 
ter M is prepared in the state \(f>) which we assume to be 
the Gaussian 46] 



da; 



(1) 



in the scaled position x with vanishing first moments. 
In the second step system and meter interact. At this 
point the interaction will be simply modeled by the von 
Neumann Hamiltonian 14 711 



H = gS(t - t )A s PM, 



(2) 



where g is the interaction strength and As represents a 
certain observable of the system which is coupled to the 
momentum operator j?m of the meter. The interaction 
happens instantaneously at time to and is weak, i.e. g is 
small. In Sec. IIVI all these model-like assumptions are 
shown to be implementable in a quantum optical setup 
for state reconstruction. 

A postselective measurement is then performed in the 
third and fourth step. Postselection begins by making a 
projective measurement of an observable B$ of the sys- 
tem. Wc obtain one of its eigenvalues bj with probability 



P(*i) = <*l (\ b j)s(bj\® 



M 



I*) 



(3) 



where \^) = exp[—igAsPM]\ip)\<fr) denotes the joint state 
of system and meter after the interaction. 

The fourth and final step consists of measuring a suit- 
able observable 6 m of the meter, given that we have 
found the eigenvalue bj before. The corresponding post- 
selected expectation value 



E w [, 



1 



om\ 



p{b 3 



(\bj)s(bj\ ® 6m) I*) (4) 



can now be approximated by keeping all terms up to 
linear order in the coupling g. With the notation [, ] for 
the commutator and {,} for the anticommutator, this 
yields H 



E 



U) 



6 M } =igRe [A^ } ] (0 1 [pm , 6 M ] 1 4>) + 
+ gTm[A$](4>\{p M ,d M }\<l>), 



(5) 



where we have used the vanishing first moments of 
Eq. (fTJ) , and the definition of the weak value [l| 



(6) 



The use of weak measurements for quantum state recon- 
struction [H, I34U361 . |45{ relies on the special form of this 
definition. Notably, it resembles expressions as they are 
used in interferometric state reconstruction algorithms 
48]. Before using this special form of the definition of 
the weak value to derive reconstruction relations in Sec. 



IIII1 we finally look at the case where the observable 6m 
to be measured in step 4 is the scaled position x m of the 
meter. Using the relations 



{<j)\{pM,X M }\4>) = 



and 



[P 



M,X M \ = -1 



(7) 



we can simplify the expectation value Eq. ([5]) and arrive 
at 



E^[x M ] =gRe 



si*,, 



(8) 



So, for this special choice of meter observables, the real 
part of the weak value Aw is a measurable quantity. 



III. RECONSTRUCTION RELATIONS 

In this section we will present the very heart of our re- 
construction scheme. It turns out that the real part of a 
weak value as in Eq. (|8]) and the postselection probability 
Eq. ^\) can be used to reconstruct an unknown state \ip) 
of the system. We will focus on the case leading to re- 
construction relations for the momentum representation 



(P|^)=^(P) = ^(P)|. 



MP) 



(9) 



Note that this means to reconstruct the modulus |^>(-P)| 
and the phase <p(P). This also means that the relations 
which wc will present are not directly applicable to the 
case of a mixed state. 

We start by regarding a weak measurement of the posi- 
tion observable As — X$ and a projective measurement 
of the conjugate variable, the momentum B$ = Ps, in 
the postselection step. If we insert these special choices 
into the definition of the weak value, Eq. ©, we get 



X (P) 



p 



(10) 



where we have also used the momentum representation of 
the position observable Xs- Finally, rewriting this with 
the polar decomposition Eq. (|9]) results in 



X (P) 



MP)\ 



MP) 
dp ' 



(11) 



An analogous relation was first mentioned in the context 
of Bohmian mechanics and its connection to weak values 
HH EIL Ell • We will utilize it to derive a reconstruction 
relation which then turns out to be particularly suited for 
the quantum optical implementation in the next section. 
By taking the real part of Eq. (|lll) we find in particular 
that 



Re 



X (P) 



MP) 
dp ' 



(12) 



If we now again use the special choice of observables As = 
Xs and P>s = Ps in Eq. © and combine the result with 
Eq. (|12l) . we arrive at 



MP) 
dp 



1 



E (p) 



X M ■ 



(13) 



This leads directly to the reconstruction relation for the 
phase, 



s-. m 



3 



<P(P) 



dP' E< p ' 



[XM 



(14) 



in which the lower integral bound is arbitrary, since any 
global phase of the state |?/>) cannot be detected. 

To complete the state reconstruction we still need the 
absolute value IVK-P)! of the system state in momentum 
representation. Hence we calculate the probability Eq. 
([3]) of finding the eigenvalue P in the projective measure- 
ment on the system. When we again retain only terms 
up to linear order in g, we arrive at 

p(P) = (tf | (\P)s(P\ ® 1m) I*) = |^0P)| 2 , (15) 

keeping in mind that the first moments of the meter state 
\4>), Eq. ([I}, vanish. Thus the reconstruction relation for 
the absolute value reads 



MP)\ = VpiP)- 



(16) 



Note that the information on the modulus IV'C-P)! just 
comes from simple projective measurements of momen- 
tum performed in step 3 of the scheme shown in Fig. [TJ 
The corresponding information on the phase ip(P), ac- 
cording to Eq. (|14|) . can be extracted only from a set of 
postselected data, i.e., a continuum of expectation values 
ordered with respect to the momentum of the system. 



IV. QUANTUM OPTICAL IMPLEMENTATION 

The reconstruction relations introduced in the previous 
section rely on the special form of the interaction Hamil- 
tonian Eq. ([2]) . So it is a crucial question to ask where in 
nature the momentum observable of one subsystem, here 
Pm of the meter M, is coupled to the position observable 
of another subsystem, here X$ of the unknown system 
S. In the next section we recall a straightforward optical 
realization for which the derived reconstruction relations 
hold. For another optical version with a very different 
aim, see Ref. 1501] . 



A. Beam splitter interaction 

In linear quantum optics [51j a beam splitter simply 
couples two modes of light; see Fig. [2j One mode rep- 
resents the system S prepared in the state \ip) and the 
second mode will be the meter in state \<p). The joint 
state after the beam splitter reads [3(| HH 



-i 6 (XsPm - Psxm) \1>)\4>), (17) 



1*) = exp 



where 9 is the parameter that determines reflection and 
transmission of the device [521 ■ Here the observables 



V2 



and Pa 



iV2 ( 



■,t 



(18) 
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FIG. 2. Schematic drawing of a possible quantum optical im- 
plementation. The system mode S and meter mode M are 
prepared in their respective states and get coupled by a beam 
splitter (BS) leaving them in the joint state Afterwards, 
using homodyne detectors HDS and HDM, momentum and 
position quadratures are measured. Based on the momentum 
quadrature readings on HDS we can postselect the expecta- 
tion values of the position quadrature at HDM. 



are the position and momentum quadrature observables 
[o"H of the system mode, which are defined by the usual 
creation and annihilation operators a s and as- Analo- 
gous relations hold for the quadratures xm and pm of 
the meter mode. Hence we see that the beam splitter 
interaction closely resembles the interaction modeled by 
the Hamiltonian of Eq. © . The additional term, deter- 
mined by the product Psxm, will actually not affect the 
postselected expectation value Eq. ([8]). 



B. Implementation 

We can now specify the calculations in Sec.[H]using the 
beam splitter interaction given by Eq. (JT7J). The state 
of the system mode is still arbitrary; the special state of 
the meter mode is given in Eq. ([T]), where |x) is now a 
position quadrature eigenstate. Thus Eq. ([I} is just the 
x representation of the vacuum state \<f>) = |0). 

The first step of the postselection is now a projective 
measurement of the momentum quadrature observable 
P>s = Ps, as depicted in Fig. [5] The probability p{P) 
of finding a certain eigenvalue P can then be calculated 
based on Eqs. (fT5j) and (fl"7|) . When we neglect all terms 
of order higher than linear in 6 and keep in mind that 
the first moments of |0) vanish, we will get the result 
p(P) = \ip(P)\ 2 in complete analogy to Eq. (|T5|) . 

Now we can continue with the expectation value of the 
position quadrature x m , also neglecting all terms of order 
higher than linear in 8, which results in 



E^[x m ] = 6{rc [XW 



— Im 



P (P)] ) 



(19) 
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Furthermore we find 



Im 



= Im 



(p\p s m 



= 0, 



(20) 



since Ps is a Hermitian operator. Thus the expectation 
value Eq. ([T9|) simplifies and is now equivalent to Eq. (j8]) 
when we replace the interaction constant g by the beam 
splitter parameter 9. The additional term Psxm in the 
beam splitter transformation Eq. (|I7[) . from which the 

(p) 

imaginary part of the weak value Pq, originates, has no 
influence on this expectation value. Hence a beam split- 
ter interaction with a small parameter and a specifically 

chosen postselection measurement can be used to obtain 

( p) 

the real part of the position quadrature weak value Xw . 
Based on this information and the knowledge about the 
probability p(P) of finding an eigenvalue P in the post- 
selection measurement, the unknown state \ip) of the sys- 
tem mode S can be reconstructed using exactly the re- 
construction relations derived in Eqs. (|I4[) and (|I6[) . 

Note that the reconstruction relations for phase, Eq. 
(|I4p . and modulus, Eq. (|I6I) . of a momentum represen- 
tation rely on the special choice of the weakly measured 
observable As = Xs and the postselection observable 
Bs = Ps as well as on the form of the interaction, Eq. 
(TT7|) . A different choice of observables and interaction 
may thus lead to reconstruction relations suitable for rep- 
resentations in terms of other degrees of freedom. This 
suggests that the general idea behind this reconstruction 
scheme might be applicable to a broader range of phys- 
ical systems, not only light modes coupled by a beam 
splitter. 



V. SIMULATION 

To see if the reconstruction scheme presented in the 
previous paragraphs is really feasible, we will numerically 
simulate the setup shown in Fig. [2] and investigate the 
various influences of relevant parameters. 



A. Influence of 6 

We will first focus on the errors introduced by neglect- 
ing higher orders of the beam splitter parameter 9. More- 
over, it is then interesting to see how a statistical error 
in the measured observables due to a finite number of 
measurement runs influences the reconstruction quality. 
To give a first example, we choose the state 



= AA(|a = I) + |/3 = 2e i7r4 / 5 )) 



(21) 



of the system mode S, where J\f is a normalization con- 
stant and | a) and \(3) are coherent states of light [53ll54j|. 
Note that this state in fact contains nontrivial phase and 
modulus dependencies which allow us to simulate essen- 
tial features of the reconstruction. 




(a) exact 




(b) reconstruction 1: 9 = 0.05, 8 = 0.0003 




(c) reconstruction 2: 9 = n/4, 8 = 0.34 

FIG. 3. (Color online) Error introduced by just keeping terms 
up to linear order in 9 in the weak- value expressions. Depicted 
are the corresponding Wigner functions: W e of the exact 
state, Eq. (|2I[) . Wi of the reconstructed state \ip rec (9 = 0.05)) 
and Wi of the reconstructed state \iprec(fi = 7r /4))- The re- 
construction error S, Eq. I)22[), quantifies the total error. 
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First we have to decide how small 9 shall be chosen 
to justify neglecting all terms of order higher than linear 
in 9. Therefore we calculate the exact expectation value, 
Eq. Q, and the exact probability, Eq. by using the 
full beam splitter transformation, Eq. (|17|) , applied to the 
system state \ip), Eq. (|21|) . and the meter state \(j>) = |0). 
Next we insert the results of these calculations into the 
reconstruction relations given in Eqs. (TT4")) and (|16p . As 
these relations have been derived disregarding all terms 
of order higher than linear in 9, they might be rough 
approximations, unless 9 is small enough. To quantify 
the error, we compute the reconstruction error 



6 = l-\(ip\ip rec (0)}\'' 



(22) 



between the reconstructed state \ip re c(9)} and the exact 
state 

In Fig. [3] we compare Wigner functions [55l - f57| of the 
reconstructed states to the Wigner function of the ex- 
act state and we also show the corresponding reconstruc- 
tion error, Eq. (|22l) . It is clearly visible that the choice 
9 = 0.05 for the beam splitter results in a small error 
introduced by retaining only terms of order up to lin- 
ear in 9. Moreover, we also note that even for a 50:50 



beam splitter [cf. Fig. 3(c) one can still reconstruct some 
features of the state 



), Eq. (|21j) . by using the proposed 
scheme. However, negative parts of the Wigner functions 
are very sensitive to a badly chosen parameter 9. 

Remarkably for a 50:50 beam splitter, i.e., 9 = tt/A, 
the setup depicted in Fig. [2] is nothing else than an eight- 
port interferometer (see, for example, [HI]). In this case 
the reconstruction error S is large for our way of process- 
ing the measured data. Yet the Husimi-Kano Q function 
p&L l59l ] may be obtained directly as the joint probabil- 
ity distribution of such an eight-port interferometer [60l ] . 
However, it is quite complicated to extract the underly- 
ing quantum state from a Q function. On the other hand, 
if 9 is small we can directly obtain the absolute value and 
the phase of an unknown state by using the simple rela- 
tions Eqs. (fl4|) and (fTBj) . At this point it is interesting to 
see how in two different measurement regimes two very 
different ways of processing measured data each lead to 
the full state information. 



B. Statistical error 

Up to now we have dealt with exact expectation values 
and full probability distributions, which in a real experi- 
ment are accessible only in the limit of an infinite number 
N of measurement runs. To account for the influence of 
the statistical error due to a finite N, we also perform 
a Monte Carlo simulation for the optical setup shown in 
Fig. [2] with the state of Eq. (|2ip. For this simulation, as 
for any real experiment, we furthermore have to choose a 
certain binning AP [6lj of the momentum quadrature 
values P. Hence we obtain approximated probability 
distributions for the P quadratures and approximated 
values for the postselected expectation values Eq. (TH)|) . 



These estimates for probability and expectation value are 
then inserted into the reconstruction relations Eqs. (| 14|) 
and (fT6j) . 

Figures [4] and [5] represent two typical measurements 
for a properly chosen beam splitter parameter 9 = 0.05. 
The number of measurement runs is N — 10 4 and 
iV = 10 5 , respectively. We show the pure measurement 
data, namely, approximated probabilities and approxi- 
mated expectation values, as well as the resulting recon- 
structed phase. The steps in the probabilities and ex- 
pectation values as well as the edges in the reconstructed 
phases are an artifact of the finite bin width AP. As 
expected, increasing the number of measurement runs N 
improves the reconstruction quality. 

Furthermore, we notice that in general those parts of 
the reconstructed phases located in the vicinity of P = 
are in better agreement with the exact phase than outer 
parts situated at larger values of P. There are two rea- 
sons for this. First, according to Eq. ([14]). we start the in- 
tegration of the postselected expectation value at P = 0. 
Therefore, the reconstructed and exact phases agree per- 
fectly at this point. Then we need more and more ex- 
pectation values to reconstruct the phase for increasing 
values of P. Hence the statistical errors involved in ev- 
ery single expectation value sum up. Second, the statis- 
tical error of a single expectation value increases when 
the corresponding probability of obtaining the outcome 
P decreases. This is most clearly visible from the outer 
parts of the estimated expectation values in Figs. S]and[Sl 
which tend to fluctuate more strongly around the exact 
value as the probability p(P) gets smaller. 

Due to these intricacies, the reconstructed phase is 
more sensitive to statistical errors than the reconstructed 
modulus. Even for N = 10 measurement runs (see 
Fig. [5]) the reconstructed phase does not match the ex- 
act value as nicely as does the reconstructed modulus. 
For N — 10 (see Fig. [4} this contrast is even more pro- 
nounced. As in other reconstruction schemes [30ll3lll48j. 
it is harder to reconstruct the phase than the absolute 
value of the quantum state. 

However, in the scheme discussed here this sensitivity 
of the phase can be fully traced back to the measured 
data: The set of postselected expectation values needs to 
be known with good statistical confidence and on the 
whole P axis. This also means that quantum states 
\ip) with regions in the quadrature distribution |^(P)| 2 
that are nonzero and separated by an intermediate region 
where the quadrature distribution is equal to zero cannot 
be reconstructed easily. In these gaps we will never find 
a postselected expectation value and hence we will lose 
the phase relation between parts of the quadrature rep- 
resentation ip(P) separated by these gaps. To a certain 
degree this can be seen already in Figs. 0] and [5j For some 
values of P we are missing the corresponding expectation 
values (see, e.g., P = —2.75 in Fig. [5J, and hence we lose 
track of the phase. However, as this happens here in a 
region where |?//(P)| stays exponentially small, it is not 
very problematic and the reconstruction error 8 is still 
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FIG. 4. (Color online) Monte Carlo simulated (red, solid) 
and exact (black, dashed) values of probability p(P), expec- 
tation value E^ p '[o;m], and phase tp(P) as a function of the 
momentum quadrature P for the state in Eq. (|21[) . Depicted 
are results for TV = 10 4 measurement runs for a beam splitter 
parameter 9 = 0.05 and a bin width AP = 0.1. The recon- 
struction error is 8 — 0.08. One can see that the statistics is 
good enough to reconstruct the modulus |i/>(P)| = ^/p(P). 
However, the postselected expectation value shows consid- 
erable deviations. This is even more pronounced in those 
regions where we have a small probability to find a certain 
value of the momentum quadrature. As a consequence the 
reconstructed phase soon differs from the exact value. 



FIG. 5. (Color online) Monte Carlo simulated (red, solid) 
and exact (black, dashed) values of probability p(P), expec- 
tation value E^'I&m], and phase <p(P) as a function of the 
momentum quadrature P for the state in Eq. (|21[) . Depicted 
are results for TV = 10 s measurement runs for a beam split- 
ter parameter 9 = 0.05 and a bin width AP = 0.1. The 
reconstruction error is S — 0.008. In comparison with Fig. 
[4] the statistics here is much better. The postselected data 
now show stronger fluctuations only in those regions where 
the probability p(P) is small. Hence the reconstructed phase 
deviates considerably only when the modulus has already be- 
come exponentially small. 
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small. 



C. Counterexample 

An explicit example of a state showing a distinct gap 
in the probability distribution |-0(P)| 2 is the state 



= jV(|a = 2i)-| J 8 = -2i)) > 



(23) 



where Af again marks the normalization constant for a 
superposition of two coherent states \a) and |/3). In 
Fig. [S]we show the corresponding Monte Carlo simulated 
data in comparison to the exact results. The probability 
|?/;(-P)| 2 around the origin is so small that there are val- 
ues of P which never appear on our meter HDS in all of 
the N = 10 5 simulated measurements. It is exactly those 
values of P where we can find no estimate for the post- 
selected expectation value. Therefore, the reconstructed 
expectation value is not continuous around the origin and 
we cannot use Eq. (TTJJ to reconstruct the phase on the 
complete P axis. 

However, as shown in Fig. |6l we can integrate the re- 
constructed expectation values in those regions where 
^(i 3 )! 2 is large enough and reconstruct a phase. To this 
end we have to choose two different starting points of 
the integration Eq. (|14p . each lying within one of these 
regions. This means we explicitly pick a certain "global" 
phase in each of the two regions. Hence we introduce an 
arbitrary relative phase shift between the two regions. 

Now let us compare this to the exact result also shown 
in Fig. [6] The exact value of the phase (p(P) jumps by 
7r at the origin and is otherwise constant. The two parts 
of the reconstructed phase show this constant behavior 
in the corresponding regions where the influence of the 
statistical error is small, i.e., where \ip(P)\ 2 is reasonably 
large. The difference between these two constant val- 
ues is exactly the relative phase shift introduced by sim- 
ply choosing a "global" phase in each region separately. 
Therefore, without any additional measurements we can- 
not determine this relative phase shift. Consequently for 
this example reconstruction fails. Fortunately, gaps in 
the probability distribution are visible from the measured 
data and therefore indicate whether we can rely on the 
reconstructed state or not. 

A possible workaround is the reconstruction of a uni- 
tarily transformed state, not showing any gaps. One ex- 
ample would be to measure suitable rotated quadratures 
on the homodyne detectors. This is equivalent to imple- 
menting a unitary transformation rotating the Wigner 
function of the input state. But this procedure will not 
dispose of gaps in a state with a rotationally invariant 
Wigner function, such as, for example, a Fock state. So 
there is no general recipe to find a unitary transformation 
that is easily implementable and removes the gaps. 
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FIG. 6. (Color online) Monte Carlo simulated (red, solid) 
and exact (black, dashed) values of probability p(P), expec- 
tation value E^'J&m], and phase <p(P) as a function of the 
momentum quadrature P for the state in Eq. (|23[) . Depicted 
are results for N = 10 s measurement runs for a beam splitter 
parameter 6 = 0.05 and a bin width AP — 0.1. The small 
probability at the origin causes gaps in the simulated expec- 
tation values. Hence we lose track of the phase at the origin 
and can only guess the relative phase between the two regions 
left and right of the origin. Therefore direct reconstruction of 
this state fails. 
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VI. CONCLUSION 

In this paper we have shown how to use weak mea- 
surements to reconstruct a continuous-variable state in 
quantum optics. The modulus of a suitably chosen rep- 
resentation can be reconstructed from projective mea- 
surements on the system itself, while the corresponding 
phase has to be extracted from a complete set of post- 
selected expectation values measured on a weakly cou- 
pled meter device. An appealing facet of the presented 
reconstruction scheme is the simplicity of the reconstruc- 
tion relations which turn measured data into quantum 



state information. Moreover, this measurement concept 
can be realized with basic elements of linear optics and 
homodyne detection. However, the sorting required for 
the postselection of data poses additional problems. The 
reconstruction of the phase is not straightforward for 
states with gaps in the probability distribution which de- 
termines the sorting. A generalization of the presented 
scheme to mixed states is desirable, but based on the 
presented reconstruction relations not straightforward. 
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